SPIN-ORBIT PENDULUM IN DIRAC OSCILLATOR * 

Marcin Turek 

O ■ Institute of Physics, Maria Curie-Sklodowska University, 

^ '■ 20-031 Lublin, Poland 

Piotr Rozmej 

Division of Mathematical Methods in Physics, Institute of Mathematics, 
^ ! Technical University of Zielona Gora, 65-246 Zielona Gora, Poland 



(N 
O 

o 



X 



Ist February 2008 



Abstract 



^ ^ The dynamics of wavepackets in a relativistic Dirac oscillator (DO) is considered. A 

^ • comparison to nonrelativistic spin-orbit pendulum effect is discussed. Particular relativis- 

tic effects, like Zitterbewegung in spin motion, are found in Dirac representation. This 
Ch i trembling motion disappears in Foldy-Wouthuysen representation. A substantial differ- 

^ ' ence between the dynamics of wavepackets corresponding to circular and linear orbits of 

a particle is obtained and discussed. 



1 Introduction 

Few years ago we were considering a motion of a wave packet (WP) representing a fermion 
in a spherical symmetric potential (e. g. the nonrelativistic t hree- dimensional harmonic oscilla- 
tor - 3dH0) with a spin orbit coupling. We have discovered a phenomenon called spin-orbit 
pendulum, in which the spin of the particle, initially polarized along a well-defined direction 
makes a reversible transition to a state in which no direction is preferred. In other words, one 
can observe a reversible transfer of the angular momentum between the spin and orbital motion. 
This behavior is very similar to revivals of population inversions within the Jaynes-Cummings 
model. 

We would like to present the relativistic version of the spin-orbit pendulum using the Dirac 
oscillator ||5|-^ instead of the 3dH0. The mechanism of spin collapses and revivals remains 
the same, but we have found several new effects caused by the relativistic regime: broken 
periodicity, the Zitterbewegung of spin mean values and the appearence of the WP components, 
that are made of negative energy states and can be interpreted within the hole theory. 

* Research partially supported by KBN grant No. 5 P03B 010 20 and Technical University of Zielona Gora 
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In contrast to the work JlSl, where the WP motion corresponding to circular orbits has 
been presented for DO, we discuss here initial conditions leading to a different topology of the 
orbital motion, i. e. hnear orbits and then compare results with that of circular ones. Results 
presented in |]13[ and in the present paper are the first ones on 3d wavepacket dynamics in DO. 
The results presented in [|T6l ,p^ correspond to the relativistic 1+1 dimensional case of Dirac 
oscillator. 



2 The Dirac Oscillator 



A system called by Moshinsky and Szczepaniak [^, [TD| the Dirac oscillator (DO) had been 
considered several times in the past in different contexts. Therefore here we only recall 
the most important relations relevant for a construction of WP and its evolution. The DO is 
described by the equation: 



ih 



'di 



Hdo^ = c[cx. ■ [p — imujr/3) + mc/3] 



(1) 



where cx.,j3 are usual Dirac matrices. It was shown that both the large and small components 
of an DO eigenstate obey equations: 
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Expression on the right side is, except an irrelevant constant, a three-dimensional harmonic 
oscillator (3dH0) Hamiltonian with a very strong spin-orbit coupling term. Hence, both the 
components must be of the form \Nljmj). The energy spectrum is given by 
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where 

2{N-j) + l forj = / + i 
2{N + j)+3 forj = /-| 

The parameter r enables, if it is small enough, a transition to the nonrelativistic limit - quan- 
tity E — mc^ gets nearly equal to the energy spectrum of the nonrelativistic 3dH0 with LS 
coupling. A requirement that \Ef; and belong to the same values of E and j results in 
the following form of normalized eigenstates corresponding to positive and negative energy 
eigenvalues, respectively 
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The following notations have been introduced: 



cjo = {rnc^)/h and tUAr^j = \E]^ij\/Ti. An exact Foldy-Wouthuysen (FW) tranformation |TT| for 
DO is also known [Q. In this representation states with negative energy disappear, which 
makes calculations much easier. The spectrum is the same as in the case of Dirac representa- 
tion m. 



3 Initial form of the WP 



We study the evolution of a Gaussian-shape WP which is initially centered at rg and has the 
average momentum po (i. e. a 3dH0 coherent state). Moreover, the initial WP is an eigenstate 
of the spin pointed at some arbitrary direction defined by a and /3, which without loosing of 
the generahty could be chosen as real numbers 



^(r,t = 0) 
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where a = ^h/muj. 

In the following we will present the behavior of two kinds of the WP 
the particular choices of initial conditions: 

• The circular WP which is obtained by choosing 

yq = xqX = (tVn X Po = Poii = cr~^VN y . 



corresponding to 



(9) 



Such WPs have already been considered in our previous papers for nonrelativistic HO [|I| 
and relativistic DO UT 



The linear WP defined by the following choise 



(10) 



Motion of linear WPs for nonrelativistic HO have been discussed in 14, 15 



Symbols x, y, z denote unit vectors along the relevant axis. The words circular and linear 
are used to emphasize different shapes of trajectories of considered coherent states. In a pure 
HO case such WPs move along circular and linear classical trajectories, respectively, without 
spreading. 
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4 Scheme of calculations 



The circular WP at t = O could be written in the basis of 3dH0 eigenstates; 



1=0 



O 

voy 



(11) 



where decomposition coeficients are given by: 



A/ 



-l)'exp(-A^/2) 



(12) 



Note, that this expression involves states with N = I = m only, taken with Poissonian weights. 
The initial form of the hnear WP is more complex 
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Introducing a complex parameter 
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we can write these coefficients as 

\ni = i-iy 



n (2/ + l)!!2'-i 

\\ {2iy.j\{N + i + i)\\ 



exp 



-KP 



C^, where j 
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Due to cyhndrical symmetry with respect to the Oz axis the WP ([TS]) contains states with 
m = O only. Next step is to obtain an exphcit form of the WP at an arbitrary time instant. It 
is convenient to perform a transformation to the coupled basis \Nljmj) because using (|^) and 
(P) we easily get a simple formula describing the evolution of the single \Nljmj) state: 
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where time-dependent coefficients A and B are defined as 
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These formulae allow us to obtain the form of the WP at given time t. However, for further 
calculations, e. g. evaluation of the spin averages, we return to the decoupled basis \Nlm)\sms) , 
which is more useful for this purpose. As a result of this procedure we obtain 



\^lin{t)) 

Every component of this bispinor has rather comphcated form 
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We have used here shortened notation: e. g. Ajv;>(t) = ANij=i+i/2{'t) ■ 
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The evolution of the WP in the FW representation is by far more simple. As previously, 
the initial form of the WP is given by (pISl). Repeating the procedure described above we get 
expressions for the WP at an arbitrary time t. Note, that the WP consists of eigenstates with 
positive energy only, so the small component - equal to zero - has been omitted for clarity 
reasons: 

i*"-w> - ( K:!: ) ■ (2^) 

Parts with spin up and down are written as 

\^i)t = Y.>^Ni{a amit) \N, 1,0) + P bmit) \N, 1,-1)} (25) 

Nl 

\(p2)t = J2>^Ni{P amit)\N,l,0) + abmit)\N,l,l)}, (26) 
m 



where 



amit) = expi-iujm>t) + expi-iujm<t) (27) 



//(/ + 1) 

bmit) = + 1 ^ expi-iujm>t) - expi-iujm<t) ) ■ (2^ 



5 Spin averages for a single |A^, /,0) [a\ t) + /3\ |)] state 



Firstly, let us discuss the behavior of the spin vector defined by the average values of spin 
operators in the case of the state 



mt = 0)) = \Nio) 
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which is one of the component states of the hnear WP ([TS]). In the FW representation we easily 
obtain: 
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where u^i = ujni> — ujnk- Hence, the vector of spin averages consists of two parts. Both of 
them he in the xOz plane. One is a constant vector which inchnation angle 9i with respect to 
the Oz axis is given by: 



tan^i 
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The second part of the spin vector oscillates with the period 2T:/ujMh The average values of 
angular momentum behave in a similar way: 

{L^)t = a(5 '^^^l^ly [ 1 - cos{uNit) ] 

{Ly)t = O (32) 

It is easy to check that the mean value of the total angular momentum (J) = (L + 1/2 cr) is 
conserved during evolution: 

{Jy)t = O (33) 
and that the length of this vector is 

\/T,^=i{Ji)t = ^|(q;2 + /32)2 ^ because for hnear WP 

L = 0. 



6 Autocorrelation function and spin averages 



An autocorrelation function at given time instant t is usualy defined as the projection of the 
WP I* (t)) onto the initial WP, that is 



A{t) = (*(0)|*(t)) 



(34) 



Hence, the autocorrelation function (or rather square of its modulus, which is a real quantity) is 
a widely used tool for illustrating recurrences in WP's behavior. As O < \A{t)\^ < 1, the value 
of |^(i)P indicates the degree of restoration of the initial shape and position of the WP. For 
our purposes we use a slightly modified (due to relativistic regime) definition of autocorrelation 
function: 

A{t) ^ I ^\f,Q)<i/{f,t)d^r. (35) 
This leads us to the foUowing expression in the Dirac representation 
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as well. The calculation of spin averages for the fuU WP in the Dirac representation is straight- 
forward but rather tedious: 
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It is easy to see that terms containing cos[(co'Ar/> — uji\fi^)t] and constant terms are dominant 
ones. 

Corresponding expressions are substantially simplified in the FW representation: 
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These expressions are just sums of (|30|) taken with weight coefficients |A 
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In Fig. 1 the evolution of average values of spin operators is presented for the circular WP 
with the spin directed initially along Ox axis (i. e. a = (3 = l/\/2). Conventionally, the circular 
orbit hes in xOy plane. For r = 0.0001 (left part) the behavior of (crj) resembles that obtained in 
the nonrelativistic description (figs. 1 and 2 of |jl|). A rapid collapse combined with a precession 
of the spin vector (defined by average values of spin operators) is observed. At time instants 
tn = {^/'^)T, where T = 27[/uj is a classical period of the 3dH0, there are spin revivals. Due to 
the quasiperiodicity of the evolution in relativistic case (nonlinear spectrum) it is hard to say if 
the spin restores with an opposite sign for n odd, as it happens in nonrelativistic case. For the 
case r = 0.0001 (weakly relativistic) the effect of Zitterbewegung is only slightly manifesting in 
{(Jz). This phenomenon is much better pronounced in the strongly relativistic case with r = 0.5 
(central picture). In this case all components oscillate more rapidly due to the interference 
among states with negative and positive energies. The Zitterbewegung is removed in the FW 
representation (right part). One can notice the fast decrease in the amplitudes of revivals for 
bigger value of r. 

Fig. 2 shows the short time evolution of spin averages for the linear WP in Dirac represen- 
tation (compare these results with nonrelativistic case in fig. 4 of 0). Note, that there is no 
precession, one observe spin coUapses and revivals only. The transition from the pure state of 
the spin in a well-defined direction n to the mixed state is especially clear in the case of the 
spin lying initially in the xOy plane {6„ = vr/2). There is a noticeable Zitterbewegung for the 
larger value of r - once again it disappears in the FW representation. The amplitude of this 
trembling motion changes with increasing angle between Oz axis and the initial spin direction, 
with the maximum for 6cr = vr/2. 

Relativistic effects, like the Zitterbewegung and an increase of the revival time with respect 
to the HO period T grow also with increasing zq and consequently with the mean energy of the 
wave packet. It is well seen in Fig. 3 that for small zq (small WP energy) the motion is almost 
periodic with time scale very close to the HO period, while for higher energies deviations from 
the HO period and imperfections of revivals become substantial. 
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7 Long time behavior 



We found the existence of fractional revivals during the long time evolution of the considered 



system. Using the prescription given by Averbukh and Perelman [|T8[ we have estimated the 
revival time (of the second order) as Trev ~ T /r for small values of r. According to Averbukh 
and Perelman's scenario the first fuU revival is expected at time (l/2)Trev preceded by a series 
of fractional revivals. Fig. 4 illustrates such revivals for the circular WP as well as for the 
linear one. For circular WP the Standard fractional revivals at t/T = 1/8,1/4,1/3,1/2 are 
well pronounced. For linear WP, due to different topology and the existence of boundary 
conditions |T^, 2D| we see spin coUapses at t /T = 1/8 and 1/4 and fractional revivals only at 



t/T = 1/3 and 1/2. 

8 Probability densities 

Initial stages of the evolution of the highly relativistic (r = 0.5) circular WP are presented in the 
Fig. 5 in spherical coordinates i) and (p. The probability density on the sphere with the radius 
corresponding to the radius of the classical orbit is displayed. The case presented in the figure 
corresponds to the initial direction of the spin parallel to Ox axis. In both representations a 
large part of the WP does not move at all because angular velocities originated from orbital and 



spin motion cancel in DO ||13[. There are also two smaller subpackets which rotate in opposite 
directions with the same absolute value of the angular velocity. These subpackets don't spread 
during their evolution and move close to the classical trajectory (|r| = xo and = 7r/2). The 
smaller one consists mainly of states with negative energy eigenvalues. This part of the WP 
disappears (the corresponding part of the probability density is shifted to the bigger one) in 
the FW representation. 

Figures 6 and 7 display the evolution of linear wave packets starting from the center of the 
coordinate system and from a turning point, respectively. The initial Gaussian WP transforms 
into spreading rings (case of the spin along O z axis, i. e. along the orbit) that move along the 
classical trajectory. The details of this type of motion have been described (for nonrelativistic 
case) in [Tl, rB,TD]. Particularly interesting is a kind of squeezing of the WP starting from 
the turning point when it reaches the center of the potential at t = (1/4)T. For other spin 
directions rings become crescents (see Fig. 8) because cylindrical symmetry is absent in the 
initial conditions. For longer evolution (not shown) one sees that the linear WP starting from 
one turning point does not reach the other. Instead, the WP is squeezed near the center of the 
coordinate system and turned back. This unusual behavior is caused by the extremely strong 
spin-orbit coupling in the DO. 

As in the previous case of the circular WP one may observe, if the value of r is big enough, two 
subpackets travelling in opposite directions. One of them is removed in the FW representation. 
This behavior is well illustrated in Fig. 9. 
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Figure 1: Time evolution of spin components average values for the circular WP with N — 20. 
Results for two different values of r (0.0001 and 0.5) are presented. 



Figure 2: Short time evolution of the spin averagcs for the linear WP in Dirac representation 
and its dcpcndnce on the initial spin direction. Two cases of r values are presented (r = 0.0001 
and r = 0.01), (cr„) is an average value of the spin projection on the initial spin direction. 



Figure 3: Comparison of the spin averages short timc evolution for different initial positions 
of the linear WP with r=0.01. Results for Dirac (Icft) and FW (right) representations are 
presented. 



Figure 4: Spin revivals in the long time evolution of the circular (top) and the linear (bottom) 
WP in a weakly relativistic case. 



Figure 5: Initial stages of the evolution of the circular WP {N = 20, r = 0.5) in Dirac and FW 
representations. Angular coordinates 'd and ip are used to present the cross-sections of 
with the sphere |r| = zo. 



Figure 6: Evolution of the linear WP starting from the center of the coordinate system 
{zo = 0). ^'^^ in the xOz planc (containing the classical trajectory) is presented. Remember 
that there is a cylindrical symmetry with respect to Oz axis. 



Figure 7: The same as in the Fig. 6. but for the WP starting from the turning point zo 
(Po = 0). 



Figure 8: Cross-sections of the hnear WP with the planes perpendicular to the classical 
trajectory. The axial symmetry is lost if the initial spin direction is not along the Oz axis. 



Figure 9: Comparison of the evolution of the strongly relativistic linear WP (r=0.5) in Dirac 
and FW representations. 
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r=0.0001 N=20 Dirac rep. r=0.5 N=20 Dirac rep. 
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